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Fractional photon-assisted tunneling is investigated both analytically and numerically for few 
interacting ultra-cold atoms in the double-wells of an optical superlattice. This can be realized 
experimentally by adding periodic shaking to an existing experimental setup [Phys. Rev. Lett. 101, 
090404 (2008)]. Photon-assisted tunneling is visible in the particle transfer between the wells of 
the individual double wells. In order to understand the physics of the photon- assisted tunneling, 
an effective model based on the rotating wave approximation is introduced. The validity of this 
effective approach is tested for wide parameter ranges which are accessible to experiments in double- 
well lattices. The effective model goes well beyond previous perturbation theory approaches and 
is useful to investigate in particular the fractional photon- assisted tunneling resonances. Analytic 
results on the level of the experimentally realizable two-particle quantum dynamics show very good 
agreement with the numerical solution of the time-dependent Schrodinger equation. Far from being 
a small effect, both the one- half-photon and the one-third-photon resonance are shown to have large 
effects on the particle transfer. 
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I. INTRODUCTION 

Recent experimental developments for ultra-cold 
atoms in optical lattices allow to count ever smaller num- 
bers of atoms. Ref. [l] uses an optical superlattice which 
is a double-well lattice [2 to multiply single double-well 
potentials with less than six atoms at each double well. 
It now is even possible to count single atoms 0. These 
developments open new possibilities to investigate few- 
particle quantum dynamics far from the mean-field ef- 
fects often encountered, e.g., in large Bose-Einstein con- 
densates (BEC). The focus of the present paper will be 
the influence of periodic shaking on the quantum dynam- 
ics of few particles in the double wells of the optical su- 
perlattice of Refs. dill]. We will use the fact that the ex- 
perimental parameters can be chosen such that tunneling 
between neighboring double wells can be discarded [5 . 
The fact that the experimental setup of Ref. [1 can real- 
ize few particles per lattice site in a very controlled way 
is based on the Mott-insulator [6 with which the exper- 
iment starts p. 

Periodic shaking is already interesting on the single- 
particle level ^ where effects ranging from destruction 
of tunneling [8] over tunneling-control [9 to population 
transfer between defects [10 and quantum scattering in 
driven single- and double-barrier systems [11 have been 
investigated. On the many-particle level, recent investi- 
gations range from transport of bound pairs in optical 
lattices [12 over many-body coherent destruction of tun- 
neling [13] and an ac-driven atomic quantum motor [14] 
to frustrated quantum antiferromagnetism [T5l. Peri- 
odically kicked systems have been investigated, e.g., in 



Refs. [16 and references therein. 

One of the interesting aspects of periodic shaking is 
photon-assisted tunneling ^3 Both integer photon- 
assisted tunneling [19 and the related effect of dy- 
namic localization [20 have been investigated experi- 
mentally for weakly interacting BECs in optical lattices. 
The experimentally observed one- and two-photon reso- 
nances [19 essentially are single-particle effects that sur- 
vive interactions. 

Fractional photon-assisted tunneling (cf. Fig. [T]) is a 
true interaction driven many-particle effect which has so 
far been shown to be a small effect both for few par- 
ticles in double wells [181 EI] and for BECs in double 
wells [22] as well as for super Bloch oscillations [23] . The 
aim of the present manuscript is to show that fractional 
photon- assisted tunneling can, in fact, be a large effect 
which can be observed with the existing experimental 
setup [1 . Strictly speaking, there are no photons ab- 
sorbed in the photon-assisted tunneling in periodically 
shaken systems: the "photons" are modulations of the 
potential on time-scales in the kilo-hertz regime. There 
are, however, resonances which correspond to integer or 
fractional multiples of the energy corresponding to the 
shaking frequency. The tunneling at those resonances 
is called photon-assisted tunneling - a half-photon reso- 
nance corresponds to one photon transferring two atoms 
to the other well, a one-third-photon resonance corre- 
sponds to one photon transferring three atoms. 

The paper is organized as follows: In Sec.[TT| the model 
used to describe the interacting particles in a periodically 
shaken double well is introduced. In order to understand 
the physics of the time-dependent system, an effective 



time-independent model is introduced in Sec. Ill Sec- 
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tion jlVj investigates the quality of the effective approach 
systematically via both analytic calculations and exact 
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FIG. 1. (Color online) Schematic drawing of fractional 
photon- assisted tunneling for the example of a 1/2-photon 
resonance for the experimentally relevant 1 case of A/" = 2 
particles. The energy of one "photon" is enough to make two 
particles tunnel. However, the "photons" are time-dependent 
potential modulations in the kilo-hertz-regime rather than 
real photons which could be absorbed. For ultra-cold atoms in 
periodically shaken double- well potentials, fractional photon- 
assisted tunneling is an interaction induced many-particle ef- 
fect. 



numeric diagonalization of the approximated model. 



II. MODEL SYSTEM 

For the description of ultracold atoms inside a double- 
well potential, a many-particle Hamiltonian in two- 
mode approximation [24 originally developed in nuclear 
physics [25j is used: 



+ h[/j.o + /ii sin(a;t)) (^4^2 ~ 4^i) • 



(1) 



The operator c^-^^ annihilates (creates) a boson in well j; 
hft is the tunneling splitting, hfiQ is the tilt between well 1 
and well 2 and hfii is the driving amplitude. The inter- 
action between a pair of particles in the same well is 
denoted by 2hf<i. Dynamics beyond this model has been 
investigated, e.g., by Refs. [26 . 

In order to describe the time-evolution of the inter- 
acting system, the Fock basis {v) = \N — v^v) is used. 
The label v = {)... N refers to a state with N — v particles 
in well 1, and v particles in well 2. The Hamiltonian ([T]) 
now is the sum of two (A/" + 1) x (A^ + l)-matrices, 



(2) 



While the non-diagonal matrix Hi is given by the 
tunneling-terms of Eq. ([T]), the diagonal matrix Hq in- 
cludes both the interaction between the particles and 
the applied potential difference. For the solution of the 
Schrodinger equation the ansatz based on the interaction 
picture [27^ 



(i/|V^(t)) = a^(t)exp 



turned out to be useful l28l. 



j {v\Ho{t)\v)dt 



(3) 



For the time-dependent amplitudes a^{t), appendix [A] 
derives differential equations which are mathematically 
equivalent to the A'^-particle Schrodinger equation gov- 
erned by the Hamiltonian ([T]): 



V2 



Here, we use the definitions o_i = and un+i = as 
well as 



and 



= -koj + - 2[iV - {2j - 1)]k , j = 1 . . . iV (5) 



(6) 



where is a Bessel function of integer order (cf. ap- 
pendix |A|. 

The Eqs. Q and ([5| include the experimentally rele- 
vant [1 case of A/" = 2, for which we introduce the abbre- 



viations CT/c = T]^^"^ and a£ = rj^^ and thus 
(j/c = —kuj + 2/io — 2/^ 



(2) 



(7) 



An important observable to detect photon-assisted 
tunneling is the time averaged transfer. 



(-^trans)^ 



NT Jo 



4 ^2 



(8) 



which is accessible to experimental measurements and 
describes the expectation value of the time averaged frac- 
tion of particles found in the upper well, if all particles 
are initially prepared in the lower one. 

In Fig. [2] a two dimensional projection of the (within 
the model ([T])) exact time- averaged particle transfer is 
shown as a function of the driving frequency uj/Q and 
the interaction hc/Q for N = 2. The dynamics display 
rich features: The solid lines mark integer and fractional 
"photon" resonances as discussed in [18 . This means the 
tilt 2hjiio is bridged by a (fractional) integer multiple m 
of the driving frequency: mhw = 2fi/io- The dashed and 
dotted lines mark resonances between the interaction Hi 
and the tilt /io plus the assistance of an integer number 
of "photons" k and i respectively, i.e. dk = (dashed) 
and d-£ = (dotted). 

Another very important approach towards resonances 
in a system governed by a Hamiltonian periodic in time 
are quasi-energies [SHl El] . They account for the prob- 
lem that for such systems neither energy nor momen- 
tum is a reasonable quantum number anymore. Ac- 
cording to the Floquet theorem, which is similar to the 
Bloch theorem known from solid state physics [32], the 
eigenstates ^(t) for a Hamiltonian H{t -\- T) = H{t) 
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FIG. 2. (Color online) Two dimensional projection of the 
(within the model ([T])) exact particle transfer ^ for the ex- 
perimentally realistic situation [l] of A/" = 2 particles in a 
double well as a function of both shaking frequency cj and 
interaction k. The transfer is averaged over TQ = 100 
with 2fii/uj = 1.8, fio — 1.5f^. All particles are initially in 
the lower well. For better visibility the fourth root of the 
transfer is plotted; bright colors correspond to large trans- 
fer. Solid lines: fractional and integer "photon" resonance, 
dashed/dotted lines: resonances between the interaction and 
the tilt plus assistance of photons. White rectangle marks 
parameter range for the quasi-energy plot in Fig. [3] The 1/2 
photon resonance is visible as a large effect in the particle 
transfer near the vertical line at cj = 6Q. 




FIG. 3. Two of the three quasi-energies for N — 2 parti- 
cles. All parameters as in Fig. |2] The inset shows these quasi- 
energies for Nk/Q = 5.3. The avoided crossing indicates the 
1/2 "photon" resonance at cj/Q ^ 6. 

have the form ^(t) = ex.p{ist / h)u{t) with the function 
u{t + T) = u{t). The eigenvalues £ of the new Hamilto- 
nian H = H{t) — ihd/dt are called quasi-energies with 
the corresponding eigenstates u{t). However, for every 
quasi-energy there can be found a whole set of further 
eigenvalues e + mhu of H. Therefore the quasi-energies 
are not unique as the choice of m G Z is free, but this 
has no influence on the physical state. 

It has been shown [30 that avoided crossings in the 
quasi-energy spectrum correspond to resonances in the 
transition probability. Thus, in the case discussed here 
avoided crossings for the right choice of parameters would 
clearly indicate resonances and support the conclusions 
already drawn from the map. In the two mode system 
with N particles there are 7V+1 quasi-energies possible in 
the range between nhw and (n + l)Jiw^ each correspond- 




FIG. 4. 1/2-photon resonance for N — 2 particles. Left: 
The nearly perfect transfer between the states occurs for the 
parameters /xo — 1.5^, uj — 6^1, k — 0.7511, 2/ii/a; — 3.113. 
The transfer was obtained by periodically shaking the double 
well. Right: The probability to occupy the intermediate state 
|1, 1) through which the transfer takes place. 

ing to a certain physical state. In Fig. |3] two of the three 
quasi-energies for = 2 have been plotted. In fact, one 
recognizes an avoided crossing for the same set of param- 
eters for which the 1/2 "photon" resonance is observed 
in Fig. (2] 

Figure [4] demonstrates that fractional photon-assisted 
tunneling, far from being the small effect of Refs. [T8l[2l]- 
|23], can in fact be a large effect: the transfer between 
both wells is nearly complete. The tunneling of exper- 
imentally realistic particle numbers oi N = 2 from the 
left to the right well occurs via occupying the intermedi- 
ate state |1, 1). While the probability is small as has to 
be expected for interacting particles [4 , it is finite and 
thus responsible for the transfer. Contrary to the case of 
STIRAP which avoids occupying unstable intermediate 
atomic states by depopulating them before populating 
them by using two laser pulses at different frequencies 
and time-dependence [33], there is no need to avoid the 
population of the intermediate level in our case. Fur- 
thermore, only a single shaking frequency is involved. 
Figure [4] also shows that the half-integer resonance is re- 
lated to the co-tunneling of the repulsively bound pairs 
of Ref. [34. 



III. EFFECTIVE MODEL 

The calculations presented in this section show, how an 
approximation similar to the rotating wave approxima- 
tion [35] leads from the set of differential equations Q to 
an eigenvalue problem, which may, for instance, be easily 
solved via numerical diagonalization, even for higher par- 
ticle numbers. For some special cases such as depicted in 
appendix [Cj the eigenvalues may also be calculated an- 
alytically for = 2. As higher particle numbers are 
treated analogously, the calculations in this section are 
initially performed ioi N = 2 particles before giving the 
general A'-particle equivalent. Furthermore, the initial 
state of the system under consideration is always |0). 

While in [18 time-dependent perturbation theory was 
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applied to an equation equivalent to Eq. Q, we here use 
a different approach: large oscillation freque ncies are dis- 
carded; this occurs in each of the sums in Eq. ( |A8| ) - which 
is the 2-particle version of Eq. Q - only the phase factor 
with the smallest frequency (7^/, and a^/ respectively, is 
maintained, according to the definition 



A = cTfe/ 



(9) 
(10) 



where k' is the integer for which 
{|— /ccj + 2/io — 2/^1 ; A: G Z} reaches its minimum 
and, analogously, I' is the integer for which 
{|— ^cj + 2/io + ; ^ G Z} reaches its minimum. 
One now has to solve the differential equations 



-Ak'e 




The ansatz 







ao{t) 
ai{t) 
a2{t) 




aie' 



-iujt 
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-i{oj — X)t 

1 

-i{uj — X — v)t 



yields a time-independent eigenvalue problem 



with 



B 



ao \ 




ao 


ai 


-( 


ai 


Ci2 J 




a2 





-Ak' 
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(12) 



(13) 



(14) 



which yields three linearly independent solutions of the 
type (12). Note that the amplitudes aj{t) are functions 



of time, while the amplitudes aj of the eigenvalue prob- 
lem are not; furthermore, as soon as three linearly in- 
dependent solutions are found this also implies that one 
knows all solutions of the linear first-order differential 



equation (11). The A/'-particle equivalent of Eqs. (12) 
and (13) can be found in appendix [bJ 

Thus, by choosing the smallest frequencies as suggested 
by the rotating wave approximation, we arrived at a 
time-independent equation which tells us via minimiz- 
ing the modulus of Eqs. (|9l and (10) and inserting the 
resulting integers into Eq. l6]) which Bessel functions will 
be relevant for the tunneling dynamics. As we show in 
the following, the effective mo d el (p!2| -(14) (and its A/'- 
particle version, see Eqs. (B1)-( [B3[ )) well describes the 
exact numerics and it is thus a tool to predict interesting 
parameter regimes for future experiments. At the same 
time, it offers a much simpler approach to understand- 
ing fractional photon-assisted tunneling than previous re- 
search [iSlET]. As expected, for vanishing interactions 



Eq. ([5| shows that only a single Bessel function remains 
thus leading to, e.g., the integer photon-assisted tunnel- 
ing investigated in Refs. [19l |22] . In general, more than 
one Bessel function will be relevant in a much simpler 
way than predicted in the approach of Ref. [18^. 



However, before testing the effective model (B1)-(B3), 
it should be noted that the choice of the smallest fre- 
quency might not be uniquely defined: rather than min- 
imizing each \r]l^^\^ it might in fact be preferable to min- 
imize the modulus of the entries on the diagonal of B, 
starting from the top. For the = 2-particle case one 
thus has: 



U Gin 



(15) 



where i" is the integer i for which 
W(}x! + ^)a; — 4/io| ;^ G Z} reaches its minimum [where 
k' is the integer defined below Eq. ([9|]. This second 
approach (to minimize the sum) will also be tested in 
this paper. 



IV. RESULTS 

The effective model allows us to understand the physics 
both via analytic calculations and by exact numeric di- 
agonalisation. We start with the condition uj/2 = 
corresponding to the case when the tilt between the in- 
dividual wells is bridged by half the photon energy. This 
implies (Jk = —uoik — 1/2) — 2i<i and = — uj{l — 1/2) ^2i<i. 

Now the condition (Jk> + = 0, which corresponds to 
the half integer resonance [18] , results in /c' + = 1 (cf. 
Eqs ([9]) and (IlO|), thus leading to the eigenvalue problem 





( Cloy 










\-\ 













-k' 



C^k' 



a2 



(16) 

Here the following definition for the tunneling frequencies 
Q has been adopted: 



Q - ^ T f^^' 
V2 \ ^ 



(17) 



Analytic solutions are derived in the appendix [C) 

In Fig. [5] both the analytic expression and the (within 
the model ([T])) exact numerical solution are displayed. 
Note that furthermore the parameters are chosen such 
that \Jo{2fii/uj)\ = I Ji(2/ii/a;)|. The approximation 
shows very good agreement with the numerics; especially 
for values close to = cj/4 the relative error, plotted in 
the inset, almost becomes zero. In experiments like [1] 
and under conditions such that the single band approxi- 
mation used in Eq. ([T]) is valid, tv/fl can experimentally 
be determined with at least 1% accuracy [36]. Thus the 
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FIG. 5. (Color online) Time dependent transfer Ptrans (|8| 
for N — 2 particles with 2iii/uj — 3.113, /io = <^/4 and t* 
denotes the time for which perfect transfer is suggested from 
the analytical calculations (cf. Eq. (C4)). We start with the 
experimentally realistic |5] initial condition that all particles 
are in the lower well (cf. Ref. [l]). Solid line: exact numerics 
within the model ([T]), dashed line: analytical approximation 
obtained by inserting Eq. (C4) into Eq. (C5). The inset dis- 
plays the relative error. 




FIG. 6. (Color online) Two dimensional projection of the 
(within the model Q) exact particle transfer for N — par- 
ticles averaged over TQ = 100 with 2iii/uj — 1.8, /xo = 1.5^. 
All particles are initially in the lower well of one of the double- 
well potentials of Ref. 1 . The 1 / 3-photon resonance is a large 
effect visible in the particle transfer near the vertical line at 



full width at half maximum obtained from our calcula- 
tions is sufficiently large to be very promising for further 
experiments. 

Figure [6] shows a map in which the (within the 
model ([T])) exact time averaged particle transfer ([s]) is dis- 
played for TV = 3 particles as a function of the coupling 
strength Nn/Vt and the driving frequency uj/VL. Again 
we see many of the features already obtained from the 
two particle solution. For two particles, a half-integer 
resonance is visible as a straight vertical line at cj = 61^ 
in Fig. |2] whereas in the three-particles plot Fig. [6] only a 
resonance moving for increasing interaction to the left of 
this frequency is observable (in agreement with Ref. [9 ). 

For the same parameters as in Fig.[6]the results for the 




FIG. 7. (Color online) Two dimensional projections of the 
time averaged particle transfer. All parameters as in Fig. [6] 
(a) Approximated transfer as given by the model ( |B2| ) for 
= 3 particles by minimizing \r]\^\ for all j = 1,2,3 sep- 
arately, (b) Approximated transfer for N — 3 particles by 
minimizing the sums | ^{^^ \ for j — 1,2,3 starting with 
low j. (c) Modulus of the difference of the approximated 
transfer displayed in (b) and the exact solution from Fig. [g] 



two approximations have been plotted in Figs, [t] (a) and 

(b) , as well as the difference between the second approx- 
imation and the exact solution, which is plotted in Fig. [7| 

(c) . 

It is cl early visible that both approximations based on 
Eq. (|B2l are fairly good reproductions of the (within the 
model Q) exact solution. However, in the case analyzed 

here the second approach, in which the sums | Yl{=i Vk'' I 
are minimized for j = 1,2,3, shows better agreement 
with the exact numerics. The probably most outstand- 



6 




FIG. 8. (Color online) Time dependent particle transfer over 
dimensionless time r for = 3, = 3, 2\i\luj — 1.8, 

cj/Q = 9.175, Nn/Q. — 7.68. solid (bl ack) line: exact numer- 
ics, solid gray (green) line: Eq. (B2) with minimizing each 
\rji]^\ separately for j — 1,2,3, dashed (blue) line: Eq. (B2) 
with minimizing the sums | \ 3 — 1? 2, 3 starting 

again from low j. 



ing feature to see this is the 1/3-" photon" resonance at 
= 911, that partially vanishes in the first approxima- 
tion and stays more visible in the second. In order to see, 
whether this feature in fact indicates that the second ap- 
proximation is for this set of parameters superior to the 
first, we further analyze the time dependent evolution of 
the system for A/" = 3. In Fig. [8] both approximations 
and the exact particle transfer are plotted. One clearly 
recognizes that the second approximation much better 
reflects the actual behavior of the system contrary to the 
first. The good agreement between the time-evolutions 
excludes accidental agreement on the level of the time- 
averaged dynamics. 

Both approximations show good agreement in the rel- 
evant parameter regime (depicted by bright colors in the 
two dimensional-projection plots) for large enough fre- 
quencies while they do not fit well for very low frequen- 
cies. This is because we deal with a high frequency ap- 
proximation, i.e. the lower the driving frequency is, the 
more relevant the discarded higher frequencies become. 

In order to decide which of the two approximations 
should be preferred for the calculation of the time av- 
eraged particle transfer we have further calculated the 
sums of the squared deviations from the exact particle 
transfer. 



(^trans)T " ( ^trans ) 



nl6 



^ (-ftrans)x' 



(18) 

across the map as a function of the driving amplitude 
2/ii/cj and for different particle numbers N = 2,5,10. 
To exclude low frequencies for which the approximations 
are both not reliable, the driving frequency ranges from 
uj/Q = 3 to uj/Q = 10 in steps of 0.025. The interaction 
strength ranges from Ni^/Q = to Nhc/Q = 16 in steps 
of 0.04. The results are displayed in Fig. [9j From this 
comparison one cannot clearly decide in favor of one of 



Q_ 

< 




FIG. 9. (Color online) Co mpar ison of the two different 
approximations based on Eq. (B2) (min: 



mmimizmg \ril^^'^\ 



separately for j — 1,2, . . . , N , sum: minimizing the sums 
I X)£=i Vk'' I loT j — 1,2, ... ,N starting again from low j) as 
a function of the driving amplitude 2/ii/a; and for different 
particle numbers which can be realized with the experimental 
setup of Ref. 1 . Displayed is the sum of the squared devia- 
tions from the exact particle transfer across the map (18). 



the two approximations. Depending on the driving am- 
plitude and particle number each of them yields better 
results for certain cases, contrary to what was presup- 
posed. 



V. CONCLUSION 

Photon-assisted tunneling was investigated for ultra- 
cold atoms in a periodically shaken optical superlattice. 
Experimentally realistic particle numbers of less than 4 
particles per double well [1 were chosen which turned out 
to be ideal to investigate fractional photon-assisted tun- 
neling; the calculations were done under the experimen- 
tally realistic condition that tunneling between neighbor- 
ing double wells can be discarded [5 . The numeric and 
analytic calculations were partially done in the effective 
model derived by combining the interaction picture with 
the rotating wave approximation. While one- or two- 
photon resonances [19 are a single particle effect, frac- 
tional photon-resonances are a true many-particle quan- 
tum effect which only occurs for interacting particles. 
Although previous investigations seem to indicate that 
fractional photon-assisted tunneling is a small higher or- 
der effect jlSl [2TH23] , we have demonstrated that for 
the particle numbers of Ref. [1 both the 1/2- and the 
1/3-resonance can be large effects visible in the transfer 
of particles between neighboring wells. For the param- 
eters leading to the complete transfer of two particles, 
the agreement between the exact numerics and the effec- 
tive model is astonishingly good. In general, the effective 
model well describes the physics of the system if the sys- 
tem is close to a resonance describing integer or fractional 
photon-assisted tunneling. 

Unlike the perturbation theory approach of Ref. [TS] . 
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the results of the present manuscript are neither re- 
stricted to short time-scales nor to fractional photon- 
assisted tunneling being a small effect. Similar to the 
rotating wave approximation, higher frequencies were ne- 
glected which will, for larger time-scales, eventually lead 
to deviations even for cases like in Fig. |5] where numeri- 
cal and analytical results agree quite well. The transfer 
via photon-assisted tunneling can be observed with an 
existing experimental setup [1]. 

Note added: Recently, a first experiment with periodi- 
cally shaken superlattice has been performed by Chen et 
at. [37] . The experiment uses the setup of Ref. [l] ; photon 
assisted tunneling is clearly visible both for single parti- 
cles and for few particles. One of the results is that the 
half-integer resonance (labeled "co-tunneling of a pair of 
atoms" by Chen et al.^ cf. Ref. [34^) is clearly visible for 
the case of two particles per double well. As predicted by 
the present manuscript, the half-integer resonance does 
indeed have a large effect on the particle transfer [37j. 
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Appendix A: Model 

A set of differential equations was derived [28] which is 
mathematically equivalent to the A/'-particle Schrodinger 



equation governed by the Hamiltonian ([T]): 

iha^{t) = {iy\Hi\iy^l)K{t)a^^i{t) 

+ (^|i7i|^-l)/i,_i(t)*a,_i(t) 



In Eq. (Al), the notation a_i(t) = aN-\-i{t) 
adopted; the phase factors are given by: 



(Al) 

, was 



K{t) = exp {i [2{N - 1 - 2u)K.t - 2fiot + 2/ii cos{ujt) /uj]) . 

(A2) 

Furthermore, there are the matrix elements of the tun- 
neling part of the Hamiltonian in the Fock basis: 



and 



(A3) 



{iy\Hi\iy-l) = -^VA^-z^ + lV^ (A4) 

To simplify the expression for subsequent integrals, we 
use the expansion in terms of Bessel functions [29^ 



^iz cos(ujt) 



(A5) 



Applying the above expansion and introducing the ab- 
breviation 



;J_ 

V2 



QJj{2fii/uj) 



(A6) 



and the frequencies 



-kuj + 2/io — 2/^ 
-iuj + 2/io + 2/^ 



(A7) 



the differential equations (Al) for a system with N = 2 



particles - the smallest particle number for which the 
interaction induced fractional photon resonances will be 
observable with the experimental setup of Ref. pLj - read: 




-Ek^kc-''''' o\ fMt) 

-EkAle^""'' -E.^^e--^' UiW 1 . (A8) 

-E^A}e'-^' J \a2{t) 



The generalization to the A/'-particle-equivalent displayed in Eq. Q is straight-forward. 



Appendix B: Effective model for N particles 



In order to derive the A/'-particle equivalent of Eqs. ( |12[ ) and ( |13[ ), we transfer the reasoning involved in the derivation 

reaches its minimum. Combining the 



of those equations to the general case: let kj be the integer for which 
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equivalent of Eq. ( 13 ) with the definition ^ leads to an eigenvalue equation, 



/ ^0 \ 

aN-i 
V Sat / 



Oi 

aN-i 
V Sat y 



with the tridiagonal Hermitian matrix 







\ 







-2 — '^Jk,e-^) 



The time-dependent solutions corresponding to eigensolutions of Eq, 









(Bl) 



(B2) 



/ «o(t) \ 
ai{t) 
a2(t) 

V CLN{t) J 



Bl[ ) with eigenvalue uji read: 
\ 



a\ exp 



a2 exp 



Y ay exp 



(-.-.if) 



■ / (1) (2)\ 



(B3) 



Appendix C: Analytic calculations 



tors (C2): 



This section demonstrates that it is possible to ob- 
tain analytic results for the particle transfer for fractional 
photon-assisted tunneling. Starting with the matrix ( 16 ) 
for the experimentally relevant [T] case of = 2 parti- 
cles, one then obtains a set of eigenvalues: 



^1=0 CJ2/3 



(CI) 



and their corresponding eigenvectors: 




1 



7(2) 



rr(3) 



(C2) 

These results are now inserted into the ansatz (Il2j and 
thus yield three linearly independent solutions v^^\t)^ 
v^^\t) and v'^^^t). Any state of the system |^(t)) can be 
expressed as a linear combination |^(t)) = biv^^\t) + 
h2V^'^\t) + h2,v^'^\t). The complex coefficients hj are 
obtained from the initial condition |^(0)) = |0). The 
time dependent transfer ([8|, can then be expressed in 
terms of the amplitudes aj{t) which depend on the 
time-dependent solutions corresponding to the eigenvec- 



2|a2(t)|' + |ai(t)p 



V 



i=i 



(C3) 



In order to achieve perfect t ransfer one finds (7^/ = 
implying a;2/3 = + expression for 

the transfer all phase factors only contain frequencies, 
which are integer multiples of CJ2/3, thus leading to a 
retrieval of the initial state for 002^ = 27r. Best transfer is 
therefore expected to occur for 



(C4) 



From the definition ([9| one finds cFk' = for /c' = and 

To account for experimentally realistic conditions, that 
is adjusting the interaction strength n/Vt via modulations 
of the depth of the wells, we investigate how sensitive 
perfect transfer is to changes in the interaction parameter 
n. It is reasonable to consider changes within the range 



9 



< K < (jj /2 as k' will remain zero. Inserting these 
values into the particle transfer ( C3 ) leads to an analytic 



expression for the transfer: 
1 



\P{0,t) - P{0J2,t) - P{0J3,t)\' 



(C5) 



PiO,t)-i^iPiL02,t) + Pius,t)) 



where 



i ex.p{—iujt) 
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^7? 



(C6) 



Equation ( C5 ) is displayed in Fig. [s] 
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